Shuttle mechanism for charge transfer in Coulomb blockade nanostructures 
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Room-temperature Coulomb blockade of charge transport 
through composite nanostructures containing organic inter- 
links has recently been observed. A pronounced charging ef- 
fect in combination with the softness of the molecular links 
implies that charge transfer gives rise to a significant deforma- 
tion of these structures. For a simple model system containing 
one nanoscale metallic cluster connected by molecular links 
to two bulk metallic electrodes we show that self-excitation 
of periodic cluster oscillations in conjunction with sequen- 
tial processes of cluster charging and decharging appears for 
a sufficiently large bias voltage. This new 'electron shuttle' 
mechanism of discrete charge transfer gives rise to a current 
through the nanostructure, which is proportional to the clus- 
ter vibration frequency. 

PACS numbers: 73.23.Hk, 72.80.Tm, 72.80.Le 

The term Coulomb blockade refers to a supression of 
the tunneling current through metallic grains embedded 
in a dielectric matrix. The origin of this phenomenon 
lies in the fundamental quantization of charge in units of 
e and it occurs in systems with small quantum charge 
fluctuations. As a result the redistribution of grain 
charges necessarily associated with a current flow can 
only be made in quantized units of e. The correspond- 
ingly quantized electrostatic charging energy, which for 
small nanoscale grains can be large compared to other 
relevant energies related to temperature and bias volt- 
age, tends to block the current. A number of different 
Coulomb blockade-based phenomena have been discov- 
ered — typically at very low temperatures — over the 
last decade or so largely due to the developments 

within nanotechnology. 

In this Letter we discuss a new type of composite meso- 
scopic structure demonstrating Coulomb blockade behav- 
ior at room temperatures j| |j| . The crucial aspect of 
these structures from the point of view of our work is 
that they contain metallic grains or molecular clusters 
with a typical size of 1-5 nm that can vibrate; their po- 
sitions arc not necessarily fixed. This is because the di- 
electric material surrounding them is elastic and consists 
of mechanically soft organic molecules. These molecular 
inter-links have elastic moduli which are typically two or 
three orders of magnitude smaller than those of ordinary 
solids M. Their ohmic resistance R is high and of order 
10 7 - If? ohm, while at the same time they are extremely 
small — a few nanometers in size. 
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FIG. 1. (a) Simple model of a soft Coulomb blockade sys- 
tem in which a metallic grain (center) is linked to two elec- 
trodes by elastically deformable organic molecular links, (b) 
Dynamic instabilities occur since in the presence of a suffi- 
ciently large bias voltage V the grain is accelerated by the 
same electrostatic force towards first one, then the other elec- 
trode. A cyclic change in direction is caused by the repeated 
"loading" of electrons near the negatively biased electrode 
and the subsequent "unloading" of the same at the positively 
biased electrode. As a result the sign of the net grain charge 
alternates leading to an oscillatory grain motion and a novel 
"electron shuttle" mechanism for charge transport. 



A large Coulomb blockade effect in combination with 
the softness of the dielectric medium implies that charge 
transfer may give rise to a significant deformation of these 
structures as they respond to the electric field associated 
with a bias voltage. We have found that self-excitation 
of mechanical grain vibrations accompanied by barrier 
deformations is possible for a sufficiently large bias volt- 
age. As we show below this gives rise to a novel 'shuttle 
mechanism' for electron transport. The point is that a 
grain (the shuttle) oscillates between two turning points. 
One of them is located near a positively- and the other 
near a negatively biased electrode. Due to the Coulomb 
blockade phenomenon an integer number of electrons are 
loaded onto the grain close to one turning point and the 
same number of electrons are unloaded close to the other 
as illustrated in Fig. 1. The result is that in each cy- 
cle the shuttle moves a discrete number of electrons from 
one electrode to the other. It follows that the current 
is proportional to the mechanical vibration frequency of 
the grain. 

One can estimate the typical frequency u of such vi- 
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brations using elastic moduli for the relevant organic 
molecules. For a metallic cluster a few nanometers in 
diameter and a molecular length of 1 nm the characteris- 
tic vibration frequency is of the order of 10 10 — 10 11 s _1 . 
Hence it can be of the order of or higher than the typ- 
ical charge fluctuation frequency v — 1/RC (here C is 
the capacitance of the metallic cluster which for room 
temperature Coulomb blockade systems is of the order 
10~ 18 — 10~ 19 F). This means that charge fluctuations 
and cluster vibrations become strongly coupled in such 
systems. 

For illustration we consider the simplest possible in- 
teresting system. It consists of one small metallic grain 
connected by elastic molecular links to two bulk leads on 
either side, as shown in Fig. la. The tunnel junctions 
between the leads and the grain are modeled by tunnel- 
ing resistances Ri(x), R,2{x) which are assumed to be 
exponential functions of the grain coordinate x. In order 
to avoid unimportant technical complications we study 
the symmetric case for which R\2 — Re ±x ^. When the 
position of the grain is fixed, the electrical potential of 
the grain and its charge q are given by current balance 
between the grain and the leads [Q. As a consequence, 
at a given bias voltage V the charge q is completely con- 
trolled by the ratio R\{x)/R2{x). In addition the bias 
voltage generates an electrostatic field £ — aV in the 
space between the leads and hence a charged grain 
will be subjected to an electrostatic force F q = aVq. 

The central point of our considerations is that the grain 
— because of the "softness" of the organic molecular links 
connecting it to the leads — may move and change its po- 
sition. The grain motion disturbs the current balance and 
as a result the grain charge will vary in time. This varia- 
tion affects the work W = aV J xq(t)dt performed on the 
grain during, say, one period of its oscillatory motion. It 
is significant that this work is nonzero and positive, i.e., 
the electrostatic force, on the average, accelerates the 
grain. The nature of this acceleration is best understood 
by considering a grain oscillating with a large amplitude 
A > A. In this case the charge fluctuations between the 
grain and the most distant electrode is exponentially sup- 
pressed when the displacement of the grain is maximal. 
As a result the grain at such a turning point gains extra 
charge from the nearby electrode, as shown in Fig. lb. 
The added charge is positive at the turning point near 
the positively charged electrode and negative at the turn- 
ing point close to the negatively charged one. The sign 
change of the charge takes place mainly in the immediate 
vicinity of the electrodes while along the major part of 
the path between them the charge on the grain is fixed 
and determined by the most recently encountered elec- 
trode. Therefore the grain acts as a shuttle that carries 
positive charge on its way from the positive to the nega- 
tive electrode and negative charge on its return trip. The 
electrostatic force is hence at all times directed along the 
line of motion and thus accelerates the grain. This state- 



ment is true also for small amplitude oscillations. The re- 
sult is an electro-mechanical instability in the sense that 
even if the grain is initially at rest a sufficiently strong 
electrostatic field will cause the grain to start vibrating. 
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FIG. 2. Charge response to a cyclic grain motion. The 
dashed arrows describe an imagined, particularity simple tra- 
jectory in the charge-position plane that allows the work done 
on the grain by the electrostatic field to be easily calculated 
and shown to be positive, hence leading to an instability: For 
times t < — r when the grain is at rest at x — —5x/2 the 
charge exchange with the positively biased electrode domi- 
nates and the grain is positively charged q = 8q; at the in- 
stant t = — r it instantaneously moves to the point x = +Sx/2 
where the charge relaxes to a new equlibrium value q — —5q; 
then the negatively charged grain instantaneously moves back 
at t — +t. During this cyclic process the electrostatic force 
acts only along the direction of the grain displacement. The 
solid line describes (realistically) the charge response at large 
oscillation amplitudes. At \x\ > L the charge exchange with 
the far lead is exponentially suppressed and the charge on the 
grain takes the fixed quantized value corresponding to ther- 
mal equilibrium with the nearest lead. As a result no net 
energy is pumped into the system. 



To be more precise one can show that for small devia- 
tions from equilibrium (x — 0, q — 0) and provided q(t) is 
defined as the linear response to the grain displacement, 
q(t) = J x(t — t')x(t')dt' , the work done on the grain 
is positive. To prove this one may consider the simple 
closed trajectory in x, q-space shown with dashed arrows 
in Fig. 2. For this trajectory it follows trivially that the 
work W is positive. Since this result can be shown ||] to 
be independent of the particular trajectory considered, 
this simplifying choice does not imply any restriction on 
its validity. The positive value of W at small vibration 
amplitudes implies that self-excitation of mechanical vi- 
brations will occur when a finite voltage bias V is applied 
to our system. 

In real systems a certain amount Q of energy is dissi- 
pated due to mechanical damping forces which always 
exist. In order to get to the self-excitation regime 
more energy must be pumped into the system from the 
electrostatic field than can be dissipated; W must ex- 
ceed Q. Since the electrostatic force increases with the 
bias voltage this condition can be fulfilled if V exceeds 
some critical value V c . If the electrostatic and damp- 
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ing forces are much smaller than the elastic restoring 
force self-excitation of vibrations with a frequency equal 
to the eigenfrequency of elastic oscillations arise. In 
this case V c can be implicitly defined by the relation 
^7 = aV c Imx(u), 

where loj is the imaginary part of the complex dynamic 
modulus. In the general case, when the charge response 
is determined by Coulomb-blockade phenomena, x is an 
increasing but rather complicated function of V and there 
is no way to solve for V c analytically. However, one can 
establish the value of V c with an uncertainty given by the 
Coulomb blockade voltage V e — e/C: 

Above the threshold voltage the oscillation amplitude 
will increase exponentially until a balance between dis- 
sipated and absorbed energy is achieved and the system 
reaches a stable self-oscillating regime |J. The ampli- 
tude A of the self-oscillations will therefore be determined 
by the criterion = Q(A). Let us now consider 

large amplitude oscillations. Then, due to the exponen- 
tially strong dependence of the tunnel resistances on the 
grain position, no net energy is absorbed by the sys- 
tem outside the region \x\ < L = In Auj/ \v < A (see 
Fig. 2). Therefore, the amount of absorbed energy in 
this limit does not depend strongly on the amplitude 
A [jl0|. This results in a saturation of the energy that 
can be pumped into the system, which is of the order 
W ~ AaVq eq L. Here q eq is the magnitude of the charge 
that obtains if the grain is in thermal equilibrium with 
one lead. In the Coulomb blockade regime, at temper- 
atures fc^T -c e 2 /C, the value of q eq is quantized in 
units of e: q eq = Ne = [VC/e + 1/2] e. To estimate the 
dissipated energy we take Q{A) — ttlujA 2 , which cor- 
responds to the simple phenomenological damping force 
Fd = —jx. This approach leads to the following expres- 
sion for the amplitude of self-oscillation when F > V c : 
A ~ XVtVC/e^uo) 1 / 2 (fi is defined in Fig. 3). One can 
see from this expression that the self-oscillation grows in 
amplitude with increasing bias voltage. 

In the fully developed self-oscillating regime the os- 
cillating grain, sequentially moving electrons from one 
lead to the other, provides a new 'shuttle mechanism' for 
charge as shown in Fig. lb. In each cycle 27V electrons 
are transferred, so the average contribution to the current 
from this shuttle mechanism is 



I = 2eNf , N = 



CV 

e 



(1) 



where / = ui /2tt is the self-oscillation frequency. This 
current does not depend on the tunneling rate v. The 
reason is that when the charge jumps to or from a lead, 
the grain is so close that the tunneling rate is large com- 
pared to the elastic vibration frequency. Hence the shut- 
tle frequency — not the tunneling rate — provides the 
'bottle neck' for this process. 

We emphasize that the current due to this shuttle 
mechanism can be substantially larger than the conven- 



tional current via a fixed grain. This is the case when 
ui 3> v. Our analysis shows that the electromechanical 
instability discussed above has dramatic consequences for 
the current- voltage characteristics of single electron tran- 
sistor configurations as shown in Fig. 3. Even for a sym- 
metric double junction, where no Coulomb staircase ap- 
pears in conventional designs, we predict that the shuttle 
mechanism for charge transport manifests itself as a cur- 
rent jump at V — V c and as a Coulomb staircase as the 
voltage is further increased. 

To support the qualitative arguments given above we 
have performed analytical and numerical analyses based 
on the simultaneous solution of Newton's equation, 



-kx — 7± + aVq, g = e\JnP„, 



(2) 



for the motion of a grain of mass m and a Master equation 
for the charge distribution function P n = Tr{5(n — n)p} 
(n is the electron number operator and p is the density 
matrix of electrons), 
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-P„ = e- x ' x T(n - l,n)P n _! + e x ' x T(n + 1, n)P n+l 
-e- x / x T(n, n + 1)P„ - e x/x T{n, n - 1)P„ . (3) 



The quantities P„ completely describe the "charge" state 
of the grain in the Coulomb blockade regime and hence 
the average current through the system 



I=£= [ T e-*/ x J2nn,n + l)P n dt 

The dimensionless tunneling rates T{n =p l,n) in Eq. 
are given by 

1 \ (CV l\ n (CV 1 
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and numerically obtained I — V curves for different 
amounts of damping are shown in Fig 3. 

A more precise (and time consuming) calculation along 
the line I sketched in Fig. 3 is shown in Fig. 4. Rather 
than solving the Master Eq. (||) we here calculate the cur- 
rent utilizing a Monte Carlo type algorithm. The non- 
monotonic behaviour of the current along this line is due 
to competition between the two charge transfer mecha- 
nisms present in the system, the ordinary tunnel current 
and the mechanically mediated current I me ch (^o > t) = 
S(x(t)—x )x(t)q(t) through some cross section at x . We 
define the shuttle current as the time average mechanical 
current through the plane located at xo — 0. This cur- 
rent together with the tunnel current for the same cross 
section is shown in Fig. 4. As the damping in the sys- 
tem is reduced the oscillation amplitude grows and the 
shuttle current is enhanced while the ordinary tunneling 
current is suppressed. In the limit of low damping this 
leads to a quantization of the total current in terms of 
2ef. 
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FIG. 3. Current due to the proposed shuttle mechanism 
through the composite Coulomb blockade system of Fig. 1. 
The current is normalized to the eigenfrequency u) of elastic 
grain vibrations and plotted as a function of normalized bias 
voltage V and inverse damping rate 7 -1 (Q = \J ae 2 /AC). 
With infinite damping no grain oscillations occur and no 
Coulomb staircase can be seen. The critical voltage V c re- 
quired for the grain to start vibrating is indicated by a line. 
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FIG. 4. Cross section along the line I in Fig. 3. The to- 
tal time averaged current consist of two parts, the shuttle 
current and the tunneling current. The time averaged shut- 
tle current is the mechanically transfered current through the 
center of the system < S(x(t))x(t)q(t) >, the remaining part 
comes from ordinary tunneling. As the inverse damping 7" 
increase the shuttle current approaches the quantized value 
In/euj = 3. The tunnel current is proportional to the fraction 
of the oscillation period spent in the middle region, |a;| < A. 
This fraction is inversely proportional to the oscillation ampli- 
tude and hence the tunnel current decreases as 7 -1 increases. 
The fine structure in the results is due to numerical noise. 

In conclusion we have performed a numerical analysis 
of the charge transport by the shuttle mechanism through 
the soft Coulomb blockade system shown in Fig. 1. The 



following results summarize the discussion above: (i) a 
dynamical instability exists above a critical bias voltage 
V c ; (ii) a limit cycle in the position-charge plane exists 
for the grain (shuttle) motion above V c leading to a novel 
'electron shuttle' mechanism for charge transport; (iii) 
The current-voltage curve has a step-like structure, a 
type of Coulomb staircase, even for a symmetric system. 
Our numerical analysis was made for a symmetric case for 
which i?i(0) = i?2(0). However all phenomena discussed 
above — with small quantitative corrections — still ap- 
pear in the asymmetric case when R\ (0) 7^ R2 (0) , which 
is of great importance for experiments. In a strongly 
asymmetric situation when i?i(0) » i?2(0) the elec- 
tromechanical instability still exist but instead of one 
threshold voltage V c a set of instability "windows" associ- 
ated with the I — V steps jll| will appear. This work has 
been supported by the Swedish TFR, KVA, and NFR. 
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